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Abstract 

In a recent paper, the canonical forms of a new multi-parameter class of Abel differential equations, 
so-called AIR, all of whose members can be mapped into Riccati equations, were shown to be related to 
the differential equations for the hypergeometric 2F1, 1F1 and 0F1 functions. In this paper, a connection 
between the AIR canonical forms and the Heun General (GHE), Confluent (CHE) and Biconfluent (BHE) 
equations is presented. This connection fixes the value of one of the Heun parameters, expresses another 
one in terms of those remaining, and provides closed form solutions in terms of p F q functions for the 
resulting GHE, CHE and BHE, respectively depending on four, three and two irreducible parameters. 
This connection also turns evident what is the relation between the Heun parameters such that the 
solutions admit Liouvillian form, and suggests a mechanism for relating linear equations with N and N-l 
singularities through the canonical forms of a non-linear equation of one order less. 



Introduction 

The Heun equation [1] is a second order linear equation of the form 

a f 1 ^ e \ , a fix — q . 

y + + — r + — )y + f iw ^ = ' (1) 
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where {a, (3, 7, 5, e, a, q} are constant with respect to x, are related by j + S + e = a + 0+ 1, and a ^ 0, a ^ 1. 
This equation has four regular singular points, at {0,1, a, 00}. Through confluence processes, equation (1), 
herein called the General Heun Equation (GHE), transforms into four other multi-parameter equations [2], 
so-called Confluent (CHE), Biconfluent (BHE), Doubleconfluent (DHE) and Triconfluent (THE). Through 
transformations of the form y — > P(x)y, these five equations can be written in normal form 1 , using the 
notation of [2], in terms of arbitrary constants {a, A, B, C, D, E, F}; for the 6-parameter GHE (1) we have 
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has two regular singularities at {0, 1} and one irregular singularity at 00. The 4-parameter BHE, 

A second order linear ODE is in normal form when the coefficient of y' is equal to zero - see the Appendix. 



v"+ ( K -x 2 + Bx + c + ^ + ^jy = o, (4) 

has one regular singularity at and one irregular singularity at oo. The 4-parameter DHE, 

„ ( A b C D A\ 

has two irregular singularities at {0, oo}. The 3-parameter THE, with one irregular singularity at oo, is 

y" + (-^x 4 + Cx 2 + Dx + E^jy = (6) 

Eq.(l), originally studied by Heun as a generalization of Gauss' hypergeometric ( p F g ) equation, as well 
as these related confluent families represented by (3-6), appear in applications in varied areas 2 . As a 
sample of recent related works, in [3] quase-normal modes of near extremal black branes are found solving a 
singular boundary value problem for (1); in [4], hyper-spherical harmonics, with applications in three-body 
systems, are developed in connection with the solutions of (1); in [5], a method of calculation of propagators 
for the case of a massive spin 3/2 field, for arbitrary space-time dimensions and mass, is developed in 
terms of the solutions of (1); in [6], parametric resonance after inflation is discussed in connection with 
the solutions of a particular form of (1). The separation of variables for the Schrodinger equation in a 
large number of problems results in Heun type equations too, typically for the radial coordinate, and also 
non-linear formulations involving Painleve type equations [7] can be derived from Heun equations regarded 
as quantum Hamiltonians. A number of traditional equations of mathematical physics, as for instance the 
Lame, spheroidal wave, and Mathieu equations, are also particular cases of Heun equations. 

The solutions for these five Heun equations are the subject of current study [8]-[15]. In this paper, 
a hitherto unknown connection between Heun equations and a single multiparameter Abel equation [16], 
known to have canonical forms solvable in terms of p F q functions [17], is shown. This connection fixes one of 
the parameters and expresses another one in terms of those remaining in the Heun equations (2), (3) and (4), 
and provides exact closed form solutions for the resulting non-trivial 4-parameter GHE, 3-parameter CHE, 
and 2-parameter BHE families. The solutions are linear combinations involving 2F1 or 1F1 functions, and 
this connection with Abel equations also turns evident what is the relation between the Heun parameters 
such that the solutions of these three families admit Liouvillian form. From these results, an alternative 
approach to finding the same solutions, by exploring non-local transformations, is derived, and some of these 
solutions are shown to match those derived in [15] using an essentially different approach. 

The multiparameter BHE, CHE and GHE equations solved in this paper are non-trivial in that they are 
irreducible, not degenerate, cases: the number of their singularities cannot be reduced and the equations 
cannot be mapped into pF^ equations through extended transformations of the form 

v^PWv ( 7 ) 

where {a, f3, 7, S, k} are constants and P(x) is an arbitrary Liouvillian function. Hence, the solutions being 
presented cannot be reduced to p F q solutions of the form 

(8) 

and it is in this extended sense that, herein, we say the equations being solved are not of p F q type 3 . 

Apart from being a way to relate Heun <-> p F q equations, leading to solutions to the former ones, this 
connection Heun ^ Abel is important in itself: Abel equations also appear frequently in applications [19]- 
[22] and through this connection it is possible to study their properties by studying those of the related 

2 For a list of applications of Hcun's equations compiled in 1995 sec p. 340 of [1]. 

3 There exist symbolic computation libraries that systematically resolve the equivalence of linear ODEs under (7) - see [18]. 
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linear equations. It is implicit in the existence of this Hcun <-> Abel link that there exists an equivalent link, 
between linear equations of order n with N and N — 1 singularities and related confluent cases, through 
single non- linear "Abel-like" equations of order n — 1. 

The paper is organized as follows. In sec. 1, some results of [16] and [17] are reviewed and a connection 
between Heun and Abel equations is made explicit. In sec. 2, 3 and 4, the restrictions that this connection 
implies on the parameters entering the BHE, CHE and GHE equations are derived, and it is shown how 
a sequence mapping Hcun — > Abel — > p F q equations can be composed to obtain transformations mapping 
Heun — > pF q equations. When the aforementioned restriction on the Heun parameters holds, these trans- 
formations lead to closed form solutions for the GHE (2), CHE (3) and BHE (4), expressed in terms of 
exponentials of integrals of 2F1 or 1F1 functions. In sec. 5, taking advantage of the results of the previous 
sections, an alternative derivation is developed leading to solutions free of integrals. A discussion around 
these results is found in sec. 6 and 7. In an Appendix, the formulas relating the normal and canonical forms 
for the Heun equations are included for completeness, as well as symbolic computation input permitting the 
verification of the solutions presented. 



1 A connection between Heun, Abel and ^F^ differential equations 

The transformations being presented, relating Heun and p F q hypergeomctric equations, were obtained by 
composing transformations which map Hcun, Riccati, Abel and p F q equations according to the sequence 

Heun — ► Riccati — > Abel — ► Abel canon i ca i — ► Riccati — ► p F q 

As we shall see in sec. 5, knowing the form of these transformations, one can re-derive them in an alternative 
way, shortcutting the step which goes through Abel equations. In this section, however, the Heun <-> Abel 
connection is kept visible: Abel equations are relevant by themselves and it was through this connection 
that the Heun <-> p F q relation being presented became evident. 
Abel equations of the second kind [20] are equations of the form 

, _ Mx) y 3 + .f2(x) y 2 + fi(x)y + f (x) 

gi{x)y + go{x) 

where the {/j, gi} are arbitrary functions and either fz{x) ^ or gi{x) 7^ 0. In [16] it is shown that, 
departing from a Riccati type equation, 

y' = h 2 (x)y 2 + h 1 (x)y + h (x), (10) 

by suitably restricting the form of the mappings hi and making use of the inverse transformation {x <-> y} 4 , 
one can construct an Abel equation, 

/ = {y - Pi){y - P2){y - pz) , n , 

(s 2 x 2 + six + s ) y + r 2 x 2 + n x + r 

where the {si,ri,pi} are constants, and s 2 =/= or r 2 7^ 0. This Abel equation is representative of a multi- 
parameter class all of whose members can be transformed into Riccati equations (10) using {x <-> y}, and 
from there into second order linear equations using the Riccati — > linear mapping [20] 

*--t£t« (i2 » 

The equations of this Abel class, named "Abel Inverse Riccati" (AIR) in [16], are then generated from (11) 
by applying to it class transformations of the form 

{x — > F(x), ^fj^gfj }, (13) 
P 2 {x) y + Q 2 (x) 



4 By {x <-» y} we mean changing variables {x = u(t), y(x) = t} followed by renaming {u = y, t = x}. 
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where {F, Pi, P 2 , Qi, Q 2 } are arbitrary mappings with F' 7^ 0, P\ Q 2 — Pi Qi 7^ 0. The relevance of the 
AIR class can be inferred from the fact that most of the Abel solvable equations found in the literature 5 are 
shown in [16] to be particular members of AIR. 

An important property of (11) is that its connection with second order linear equations, that is, its 
"Inverse Riccati" character, is invariant under Mobius (linear fractional) changes of x and y. This property 
is used in [17] to accomplish a full classification of (11) in terms of six canonical forms. For that purpose, 
through Mobius changes of y, (11) is first transformed into 

y> = PM (14) 

(s 2 x 2 + si x + s ) y + r 2 x 2 + n x + r 

for some new constants {si r^}, with P(y) equal to y (y — 1), y or 1, respectively according to whether in (11) 
there are three, two or only one distinct roots pi. As shown in [17], each of these three cases splits further 
into two subcases, and the six resulting canonical forms are solvable in terms of 2 Fi, 1F1 and 0F1 functions; 
in this way, closed form p F q solutions can be constructed for the whole AIR class. 

The key observation now is that the AIR equation (14) is also connected in a surprisingly simple manner 
to the Heun family of equations. As we shall see, by applying to (14) the transformation {x <-> y}, one 
obtains a Riccati equation, and by transforming it further into a second order linear equation using (12), 
one directly obtains the GHE, CHE or BHE Heun equations (with some restrictions on the parameters), 
respectively according to the three possible values of P(y). Since the AIR (14) admits solutions expressible 
using p Y q functions for the three possible values of P(y), the GHE, CHE and BHE Heun families which can 
respectively be derived from (14) also admit closed form solutions expressible in terms of these p F q functions. 



2 Closed form solutions for a subfamily of the BHE 

Considering first the simplest case, where the three roots pi in (11) are equal, in (14) we have P(y) — 1 and 
so the AIR equation becomes 

y> = \ (15) 

(s 2 x 2 + si x + s ) y + r 2 x 2 + n x + r 

Recalling that either s 2 7^ or r 2 ^ 0, changing variables using {x <-> y} we obtain the Riccati form 

y' = (s 2 x + r 2 ) y 2 + (s 1 x + n) y + sox + r (16) 
Using (12), this equation is transformed into the second order linear equation 

(s 2 s 1 x 2 + (s 2 r 1 +sir 2 )x + s 2 +r 2 r 1 ) „ 

V = ; y - (s 2 s x + (s 2 ro + s Q r 2 ) x + r 2 r ) y (17) 

s 2 x + r 2 

This equation has one regular singularity at —r 2 /s 2 , one irregular singularity at 00, and, by rewriting it in 
normal form, it is straightforward to verify that it is the BHE equation (4) with one of its four parameters 
fixed and two other ones interrelated. For that purpose, we note first that the case s 2 = presents no 
interest since it directly simplifies (17) to a p F q equation. Assuming s 2 7^ in (14), we take s 2 = 1 without 
loss of generality. To have the regular singularity of (17) located at 0, it suffices to take r 2 = 0, and, 
taking so = s\/A — 1, the coefficient of x 2 in the normal form of (17) will be as in (4). In summary, using 
y — > v / xcxp((x (si x + 2 n) /4)y to rewrite (17) in normal form at {s 2 = l,r 2 = 0, s = s^/4 — 1}, the 
equation becomes 



y" = I x 2 



fs\r\ \ ri 2 n 3 \ . „. 

(— - r °) X+ ^ + 2-x + ^) y (18) 

which is the BHE (4) at {B = r a - Sin/2, D = -n/2, C = -D 2 , E = -3/4}. 



3 For a collection of these see [23]. 
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The relevance of this result is in that, on the one hand, (18) is a non-trivial 2-parameter form of the 
BHE for which solutions are not known in general; on the other hand, as shown in [17], the Abel equation 
(15), from which (18) is derived, is solvable in terms of 1F1 and oFi (Kummer and Bessel) hypergeometric 
functions. Therefore, a closed form solution for the BHE (18) can also be expressed using these p F q functions. 

It is interesting, when possible, to separate the Liouvillian from the Non-Liouvillian solutions of (18), since 
the former ones represent the "special cases" , where the solution is representable in terms of known functions. 
For that purpose, the parameters {r ,ri} are redefined in terms of new parameters {a, t} according to 6 

r = -2a + sir, ri=2r (19) 
With this choice, the BHE (18) also becomes an equation explicitly depending on only two parameters {cr, r}, 

y "=^ + 2ax + r 2 + ^ + ^jy, (20) 

and (15) at {s 2 = l,r 2 = 0, s = sf/4 — 1} becomes 

y i = \ (21) 

(x 2 + s 1 x+(s 1 + 2)(s 1 -2)/4)y + 2TX-2c7 + s 1 T V ; 

2.1 Liouvillian solutions for the BHE (20) when a = ±r 

As explained in sec. 2.2, when a = ±r, the BHE (20) can be obtained from an equation of the form 
y" + J(x) y' = 0, "missing y", through a Liouvillian transformation, and so it admits Liouvillian solutions, 
computable using the relatively new Kovacic algorithm [24]. Concretely, for a = r, (20) becomes 

y"=(* 2 +n + r(2x + ^ + ^y, (22) 
the aforementioned equation "missing y" is 

2x(x + r) + l 

x 

and the transformation mapping this equation into (22) is 

y^ V^c^ x+2T Wy (24) 

Hence a general solution for (22) is 7 

e -x(x+2r)/2 ^ c *(* + 2T)/2_^ TCrfi(x + r)c -^+2r)/2-T 2 

V = 7= Ci H -= C 2 (25) 

Vx Vx 

It is important to note that computing the transformation (24), which maps an equation in normal 
form, like (22), into one that is "missing y", like (23), is not a trivial operation, and is entirely equivalent 
to computing the solution for (22). These mappings are formally performed with the aid of infinitesimal 
symmetry generators, and in the case of linear ODEs, the computation of these infinitesimals indeed requires 
solving the ODE itself [26]. Besides the power of Kovacic's algorithm, which computes these solutions 
systematically just by assuming the solution field (Liouvillian), it is also remarkable that the condition 
cr = ±t for the existence of these Liouvillian solutions of the BHE (20) is directly evident in the canonical 
form of the corresponding Abel equation (29) below. 

Finally, when in (20), a = — r, the same treatment with Kovacic's algorithm results in the general solution 

e x(x-2r)/2 ^ e -x(x-2r)/2_ 1TCrf(l _ T)eI Mr)/2 + r 2 

y = -j= Ci H -= C 2 (26) 

Vx Vx 



6 Thc motivation for this particular choice of {a, r} becomes clear below, in connection with the form of equation (29). 
7 In (25) and (26), erf and erfi are respectively the error and imaginary error functions - sec [25]. 
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2.2 A solution in terms of 1F1 functions for the BHE (20) when a 2 ^ r 2 

A transformation relating (20) to a 1F1 differential equation, providing a solution for (20) when a 2 ^ t 2 , is 
constructed by composing three transformations: one which maps the BHE (20) into the AIR equation (21); 
one which maps (21) into an AIR equation admitting 1F1 solutions; and finally, one which maps that AIR 
equation into a 1F1 equation. 

Reversing the transformations used to derive (18) from (15), the transformation mapping the Heun 
equation (20) into the Abel AIR (21) is 

(J xyy' dx+^y 2 +r y) | 

T , ) < 27 » 

According to [17], the transformation mapping (21) into a canonical form of AIR admitting a p F q solution is 

-l 

2 

The resulting AIR canonical form is 




£-1, V^-°) (28) 



V ' xy + x 2 + (a 2 -T 2 )/2 (29) 

This form turns evident the motivation for introducing {a, r} according to (19). For a = ±r, the independent 
term in the denominator of (29) cancels, and hence, when transforming this equation into a second order 
linear equation, we will obtain one of the form y" + J(x) y' = with rational J{x), admitting a constant 
for solution. Since the transformation of such an equation into a normal form like (20) is Liouvillian, the 
normal form of the equation will admit Liouvillian solutions. 

After having determined a condition for the existence of Liouvillian solutions, for the purpose of relating 
the Heun equation (20) to a p F q equation, a simpler derivation is possible if instead of using (28) we use 

(30) 

which docs not lead to the canonical form (29), but still leads to an AIR equation admitting p F q solutions 



y' = — , —_ , _, 2 — (3i) 



i 

(2 x - 1) y + 2 (~ + r) x 2 - 2 r x 
Following [17], this equation is transformed into a Riccati equation, then into a linear second order one using 



leading to 



'■»- -wUi > (32) 



y" = 2(x-T)y' + 2(T + <j)xy (33) 

Finally, using {x — » (x + a) 2 , y — > {{x + a) c~ x( - T+a ^)~ 1 y}, equation (33) is obtained from the confluent 1F1 
hypergeometric equation 

xy" + (u-x)y' - fiy = (34) 

at {fi = (2 + t 2 — cr 2 )/4, v = 3/2}, from where the solution to (33), in terms of the Kummer M and U 
functions 8 [25], is 

8 An equivalent form of this solution in terms of 1F1 functions is 

V = c-^ T + CT ) (1F1 ((r 2 - a 2 )/4; 1/2; (x + a) 2 ) Ci + (x + a) 1F1 (l/2 + (r 2 - a 2 )/4; 3/2; (x + a) 2 ) C 2 ) 
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V = e"^ +CT) (M \, (x + a) 2 ) Cx + U ±, (x + a) 2 ) C 2 ) (35) 

Summarizing, we depart from the Heun Biconflucnt equation (20) and arrive at the 1F1 equation (33) 
with solution (35) through a process of the form 

Heun — > Abel -» Abel lFlsolvab i e -> ^ 

The three transformations used, (27), (30) and (32), can be combined into a single transformation, mapping 
the BHE (20) into the 1F1 (33) in one step: 

^ cx p(y^ tx + 2{t + (t) Jv dx ) (36) 

Therefore, a closed form solution for the BHE (20) when a 2 ^ r 2 is given by this transformation (36), where 
in the "right-hand-side" the value of y is given by (35). By exploring some properties of linear differential 
equations discussed in sec. 5, it is possible to express the solution (36) as a linear combination of 1F1 functions 
with non-constant coefficients and entirely free of integrals - see (85). 

An independent verification of (36), or its form free of integrals (85), as well as of the Liouvillian solutions 
(25) and (26), was performed in the Maple symbolic computation system - see the Appendix. 



3 Closed form solutions for a subfamily of the CHE 

A confluent family of Heun equations can be derived from (11) when, among the three roots Pi, only two 
are different. Hence, in (14) we have P(y) = y and the starting AIR equation is 



(s2 x 2 + si x + s ) y + r 2 x 2 + n x + r a 

As in the previous section, changing variables using {x <-» y}to obtain a Riccati type equation, then using 
(12), we obtain the linear equation 



(s 2 six 2 + (sir 2 + s 2 r 1 )x - r 2 + r 2 n) (s 2 s x 2 + (s 2 r + s r 2 ) x + r 2 r ) 

y = 1 ; ; — v y - 1 2 y 38 ) 

x{s 2 x + r 2 ) x z 

This equation has two regular singularities at {0, — r 2 /s 2 } and one irregular singularity at oo, and by rewriting 
it in normal form, its confluent Heun type becomes evident. As shown below, the implicit restrictions in (38) 
as compared to the most general case (3) consist of having one of the five parameters, E, fixed and another 
one, B, being a function of the remaining three. 

To derive the relation between the parameters of the CHE equations (3) and (38), and then a solution 
for (38) in the non trivial cases, we start by noting that, when s 2 = 0, (38) simplifies to a p F q equation; 
the interesting case is s 2 ^ 0, which is equivalent to taking s 2 = 1 in (37). The regular singularities of 
(38) are fixed at {0, 1} by taking r 2 = — 1, and the term independent of x in the normal form of (38) 
is fixed to be A, as in (3), by taking sq — s\/A + A. In summary, rewriting (17) in normal form, using 
y — ► a;( ri_1 )/ 2 sfx — lexp(si x/2) y, at {s 2 = 1, r 2 = — 1, s = s 2 /4 + A}, the equation becomes 



s 1 2 /2 + s 1 (r 1 -l)-r 1 -2r + 2A+l | si + n - 1 n 2 + 4r - 1 3 \ 

~ A+ 2~x + 2(x-l) + Ix~ 2 + T{x~lf V m 



which is the CHE (3) at 



B = -A-D-C 2 , C= l Sl 2 T \ D= 1 -^-r , E = -\ (40) 



7 



As shown in [17], the Abel AIR equation (37) can always be solved in terms of 1F1 hypergeometric 
functions, from where a closed form solution expressed using 1F1 can also be constructed for the CHE 
(39). In order to separate the Liouvillian special cases of the solutions of (39) from the generally non- 
Liouvillian ones, the parameters {r , r{\ are redefined in terms of new parameters {cr, r} according to {r = 
(1 — 2a) p 2 — si 2 /4 + sirp, r\ = 2rp — s\}. This redefinition is derived as in the previous section, from the 
canonical form of the AIR (42) (see comments after (29) and also [16]). Introducing also A = — A 2 (see [1] 
p. 94), (39) becomes 



! /-H^ + 2( "- 1)A2 - TA + 1/2 + ^ + (r2 - 2 " + 2 " A2 - 1/4 + ' ,), (41) 
x x-l x 2 4 (x - 1) ' 

that is, a 3-parameter equation, and the AIR equation (37) becomes 



V (x 2 + s 1 x-\ 2 + s 1 2 /4)y-x 2 + (2T\-s 1 )x+(l-2a)\ 2 + s 1 T\-s 1 2 /4: ^ 
3.1 Liouvillian solutions for the CHE (41) when a = ±r 

As in the previous section, we know, by construction, that at a = ±r, the CHE (41) admits Liouvillian 
solutions. Using Kovacic's algorithm, for a = r, 

(1-t)A+1/2 -Ax 

y = == (Ci + (r(2 (r — 1) A + 1, -2 A x) + 2 AT(2 (r - 1) A, -2Xx))C 2 ) (43) 

y/x - 1 

where T (of two arguments - see (6.5.3) in [25]) is the incomplete gamma function. For a = — t, 

-(1+t)A+1/2 Ai 

y = == (Ci + (2Ar(2(T + l)A, 2 Ax) - T (2 (r + 1) A + 1, 2Ax))C 2 ) (44) 

\Jx — 1 



3.2 A solution in terms of ]Fi functions for the CHE (41) when o 2 ^ r 2 

As in the previous section, a solution to (41) when a 2 ^ t 2 is constructed by composing three transformations: 
the one which maps the CHE (41) into the AIR (42); one which maps (42) into an AIR equation admitting 
1F1 solutions; finally, one which maps that AIR equation into a 1F1 equation. 

Reversing the transformations used to derive (39) from (37), the transformation mapping the Hcun 
equation (39) into the Abel AIR (42) is 



exp (- / «' + '^« 2 - < 2 * - ^XX 1 " TA + ~ + 1),2)V : <fa ) } (45) 

According to [17], it is possible to construct a transformation mapping (42) into a canonical form of AIR 
admitting a p F g solution. However, as in the BHE case, simpler expressions result if we transform (42) into 
a non-canonical AIR equation. The transformation used for this purpose is 

*4-f-A, (46) 

which maps (42) into the Abel equation 

v' = (47) 

y (2\x-l)y + 2\ 2 (T + a)x 2 -2(r + l)Aar + l V ' 

Following [17], this AIR equation can be transformed into a p F g one by combining the {x <-» y} transformation 
with transformation (12) mapping a Riccati into a second order linear equation; the combination results in 
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leading to 

., 2(x-t -1)A-1 , 2 (r + a) A 2 (x - 1) 

y" = — ?/ + — '-r* y 49 

X X 1 

Finally, using {x -» 2 Ax, i/ -» x( 1+T - vl ' 2cr+r ') A 7/}, equatio n (49) is obtained from the confluent 1F1 
hypergeometric equation (34) at {fj, = (cr — 1 + \/i — 2 cr + r 2 ) A, ^ = 1 + 2 \/l — 2 cr + r 2 A}, from where 
the solution to (49), in terms of the Whittaker M and W functions 9 [25], is 



e \x 



,(r+l)A+l/2 



(50) 



Q + (1-ct)A, v 7 ! - 2 a + r 2 A, 2 Aij d + W Q + (1 - cr) A, Vl-2cr + r 2 A, 2Ax^) Cj 

Summarizing, departing from the Heun confluent equation (41) we have arrived at the 1F1 equation 
(49) with solution (50) through a process of the form Heun — > A6eZ — » Abel 1 F lSO ivabie — ► iPi- The three 
transformations (45), (46) and (48) can be combined into one transformation, 



x 



-(r+l)A+l/2 / f(x-l)v \ 

exp Ax + 2A 2 (r + a) / V / y rfx , (51) 

which maps the CHE (41) into the 1F1 equation (49) in one step. A closed form solution for the CHE 
(41) when a 2 ^ r 2 is then given by (51), where on the "right-hand-side" the value of y is given by (50). 
Like (36) in the BHE case, this solution (51) can also be expressed as a linear combination of 1F1 functions 
with non-constant coefficients and free of integrals - see (89). An independent check for correctness of these 
solutions was also performed using symbolic computation software. 

4 Closed form solutions for a subfamily of the GHE 

Solutions for the GHE (3) are obtained from (14) by taking P(y) = y (y — 1), that is, departing from 



yjy - 1) 

(s 2 x 2 + si x + s ) y + r 2 x 2 + n x + r 
The steps to construct these solutions are the same as those of the previous sections. Using 



«' - ~. ^T— " ~ (52) 



x (1 — x) y' , ,„„. 

/ , \ (53) 

(s 2 x + r 2 )y 



the AIR equation (52) is transformed into the linear ODE 



„ _ (s2 (si - 1) x 2 + ((si - 2) r 2 + s 2 r\) x + r 2 (1 + n)) , (s x + r ) (s 2 x + r 2 ) 
J ~ x(s 2 x + r 2 ) (x- 1) V x 2 (x- l) 2 



This is a Heun equation of the form (1), with four regular singularities at {0, 1, — r 2 /s 2 , oo}. For s 2 = 0, 
(54) simplifies to a Gauss equation with 2 Fi solutions. When s 2 ^ 0, in (52) one can take s 2 — 1, and, 
putting r 2 = —a, the singularities of (54) are fixed at {0, l,a, oo}, resulting in a non-trivial Heun family 
depending on four parameters. The relation between the Abel parameters {r^, Si} and the Heun parameters 
{A, B, D, E, F} is obtained by rewriting (54) in normal form and comparing coefficients: 

9 M(a,M) = 2 i, + 1 /2 e -z/2 1 F 1 (l/2 - a + b; 1 + 26; z). 
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A = r - (2r + s )a- ^ ( (n + 1) [si + - ) + n 2 - 1 ) , 



2 V V a 

r 1 +s 1 o + l- 2( J 4a 2 + (l-i)a) 
2 a (a - 1) 



1 -r x 2 

D = — r a, (55) 

£ = (l-a)f(A + B) 2 a-(A + B)(A + B-l) + ^-^-^ 



a a" 



3 

"I 



So, with respect to the most general case (2), the restriction in the GHE (54) under consideration consists of 
fixing one of the six Heun parameters, F, and expressing another one, E, as a function of those remaining. 
To derive a solution to (54), the equation is first written in normal form and the parameters {so,ro,ri} are 
redefined in terms of new parameters {A, a, t} using 

2 / 2 \ 

s = —. A 2 , r = ( A 2 - ^- ] a - 2cr A + sir, r 1 = 2r-as 1 (56) 



4 V 4 , 

aiming at showing that there are only four independent parameters and at separating non-Liouvillian from 
Liouvillian special cases of the solution, which will happen at a = ±r. The dependence of (56) on {r, a} is 
derived as in the previous sections, from the canonical form of the AIR (57) (see comments after (29) and 
also [16]). The new parameter A is related to the Heun parameter A through so and the first equation in 
(55), and is introduced here to avoid square roots in the transformation formulas 10 . With this notation, the 
AIR (52) becomes 

„/ = v(y-V (57) 

y (.T 2 + six + si 2 /4-A 2 )y-ax 2 + (2r-asi)x-(si 2 /4-A 2 )a-2CTA + s 1 T v ' 

and the Heun equation (54) in normal form, at {s 2 = l,r 2 = —a}, appears directly expressed in terms of 
the four irreducible parameters {a, A, a, r} as 



„ _ / 2a 2 (a- 1) A 2 - 2cra(2a- 1) A+ (2r 2 - l/2)a + r + 1/2 
^ \ ax 



2 (a (a - l) 2 A 2 - o (2 a - 1) (a - 1) A + (r - 1/2) ((r + 1/2) a - r)) 



T-a+ 1/2 



(a — 1) (x — 1) a (a — 1) (x — a) 

a 2 A 2 -2aaA + T 2 - 1/4 (a-l) 2 A 2 -2cr (a - 1) A + r 2 - 1/4 3 \ 

+ ~ 2 + (x~lf + A(x-af) V 

4.1 Liouvillian solutions for the GHE (58) when a = ±r 

By construction, as in the previous sections, for a = ±t, the GHE (58) admits Liouvillian solutions com- 
putable using Kovacic's method; for a = r, 

10 The use of A 2 = —A in sec. 3 brings the same advantage. 
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x T-aA+l/2 ( _ 1 > ) (a-l)A-r + l/2 

y = (59) 

(Ci + (B x (1 + 2 (a A - r) , 2 ((1 - a) A + r)) - a B x (2 (a A — r) , 2 ((1 — a) A + r))) C 2 ) 

where Bj; (one index and two arguments - see (6.6.1) in [25]) is the incomplete beta function. The Liouvillian 
form of the solution of (58) at <r — — t is 

r+oA+l/2 / _ 1 n(1-o)A-t+1/2 

y = - y —=L (60) 

Va — x 

(Ci + C 2 (B x (1 — 2 (a A + r) , 2 ((a -1)A + t))-oB x (-2 (a A + r) , 2 ((a - 1) A + r)))) 

The existence of solutions in terms of the incomplete Beta functions is also pointed out in [14], where a 
solution for the GHE as a power series is constructed, and an approximate solution involving a combination 
of two incomplete Beta functions is derived - see also [10]. 

4.2 A solution in terms of 2 Fi functions for the GHE (58) when a 2 ^ r 2 

As in the case of the BHE and CHE equations, a solution for the GHE (58) when a 1 ^ r 2 is constructed by 
composing a transformation mapping (58) into the Abel AIR (57) with a transformation mapping (57) into 
a second order linear equation admitting hypergeometric solutions, in this case of 2 Fi type. 

The derivation of results below follows the same path shown in the previous sections. Summarizing, the 
transformation mapping the GHE (58) being solved into the AIR (57) is 

( f ((a:+(s 1 -l)/2)y 2 + ((l-2x- 31 )a + T)j / +(( 31 /2 + x)a-T-l/2)a)y' ~| 

<X^y, y^cJ vtv-D^-v) dx \ (61) 



y^x) (62) 



The transformation mapping the AIR (57) into a p F q second order linear equation is 

A (r + a) y _ si 
x(x-l)y' 2 
and the resulting equation, of 2 Fi type, is 

(2(A-l)a; + l-2( a A + r)) 2 (x a) A (r + a) 
V x(x-l) V + xHx -lf V (63) 

Combining the transformations (61) and (62), a transformation mapping (58) into (63) in one step is 



aA+r + l/2 / 1 N(l-«)A-r + l/2 , f , _ x 

y ^x_ (xjl cxp 2A(r + a) / {X f \ - c 

\/a- x V J (x-l) 2 x 2 y' 

Hence, the solution to the Heun equation (58) to which this section is dedicated is given by the expression 
above, where, in the "right-hand-side", the value of y is given by the solution to (63), that is, 



y = x aA+T - T (x-l) (1 - a)A+E - T (65) 
( 2 Fi(S + A - T, S - A + 1 - T; 1 - 2T; x) Ci + x 2T 2 Fi(S + A + T, S - A + 1 + T; 1 + 2 T; x) C 2 ) 
where, to make the structure of this solution visible, instead of {a, r} we are using 
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£ = yj(a - i f A 2 - 2 (a - 1) a A + t 2 , T = \/a 2 A 2 -2aaA + T 2 (66) 
This solution (65) in turn is computed noting that (63) is obtained by changing variables 

y - .x T - aA - (.x - y (67) 

in Gauss' 2 Fi equation 

(x 2 - x) y" + ((/i + J/ + 1) x - p) y' + fj, v y = (68) 

taken at {/j = £ + A — T, z/ = £ — A + l — T, p = 1 — 2 T}. Hence, the same transformation (67) maps the 
solution of Gauss' equation into (65). For a solution equivalent to (64), free of integrals, expressed as a linear 
combination of 2F1 functions with non-constant coefficients, see (92). These and the Liouvillian solutions 
presented for the GHE (58) were also verified for correctness using symbolic computation software. 



5 Alternative derivation of solutions free of integrals 

In the previous sections, non-Liouvillian solutions, as well as their Liouvillian special cases and the relation- 
ship between the Heun parameters for their existence, were derived for the BHE, CHE and GHE equations 
(20), (41) and (58). The non-Liouvillian solutions (36), (51) and (64), however, have the drawback of con- 
taining non-trivial uncomputcd integrals. In this section, from the knowledge of the form of these solutions 
and exploring non-local transformations, equivalent solutions free of integrals are derived. 
We start by recalling that second order linear equations 



y" = c iy ' + c a y (69) 

where a = C{(x), can always be mapped into Riccati equations (10) back and forth. The transformation 
mapping (69) into a Riccati equation is of the form 



y 



j-J G(x)ydx 



where G{x) is an arbitrary function, and the transformation mapping a Riccati equation (10) into a linear 
equation is given by (12). It is also known that any two Riccati equations can be mapped between themselves 
through Mobius transformations of the dependent variable y with variable coefficients /, = fi(x), 

h y + h 

where /1/4 — 73/2 7^ 0. If we now transform (69) into a Riccati equation using (70), then apply the Mobius 
transformation (71), and to the resulting equation we apply transformation (12), we obtain another second 
order linear equation. The composition of these three transformations is the non-local transformation 

y -» exp - / - dx , (72) 



f 3 Qy' + f 4 ry 



where 



- /1/4 - hh £ 0, T = Co H + {fi - ci/i) h - ih' + h) fi + 0, (73) 

and, in fact, this transformation suffices to generate the whole class of linear equations from any given one. 
For example, applying (72) at /1 = fi = 1 to y" = 0, we obtain an equation as general as (69). 
The particular case of (72) at /1 = fi = and / 2 = / 3 = 1, 



It- 



y-exp r-^dx) (74) 
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is relevant to the results of the previous sections: the transformations (36), (51) and (64), respectively 
mapping the BHE, CHE and GHE into p F q equations, are in fact compositions of transformations of the 
form (74) with transformations y — > P(x) y. The linear equation obtained by applying (74) to (69) is 

y"=(^-c?)y' + c y (75) 

By applying to this equation the same transformation (74), we reobtain 11 (69). So, y in the "left-hand-side" 
of (74) represents the solution to (75) or (69), respectively written in terms of the solution to (69) or (75), 
represented by y in the "right- hand-side" of (74). As we shall see in the following subsections, this is 
indeed the mechanism by which solutions to Heun equations (here represented by (75)) were expressed as 
exponentials of integrals of solutions to p F q equations (here represented by (69)) in the previous sections. 

With this understanding of matters, however, it is possible to show that, despite the integral sign entering 
(74), the solution to (75) can be derived without performing any integration. For that purpose, we note 
that, given a generic linear ODE in y, it is always possible to construct the generic linear ODE of the same 
order satisfied by p = y' . Concretely, if y satisfies (69), the ODE for p is 12 



P ' ' = {^; +Cl ) P ' + ( C1 ' + C °- CJ ^ L ) P (76) 
and by substituting y = Jp dx into (69), we obtain 

y=[ P dx = P ^ P (77) 

Therefore, when (76) can be solved, the solution y for (69) can be obtained directly from p by differentiation. 
For instance, writing p in terms of some / = f(x) and g = g(x), as 

p = fC\+gC 2 , (78) 
a solution for (69) is computed from p, without using integration, as 

y= f —^C 1 + 9 —^C 2 (79) 
co c 

Composing now the non-local transformation (74) with the introduction of p = y' , that is, plugging the 
coefficients of (75) into (76), we obtain 

p = \—-c 1 )p+\*-c l + _ 2(-j jp, (80) 

and the key observation is that the normal form 13 of (80) is the same as that of (69). Consequently, if (69) 
is of p F q type 14 , then (80) is too, even when (75) may not be (in what follows it will be of Heun type), and 
the solution to (75) can be expressed not just using the integral form (74), but also using (79) as a linear 
combination, with variable coefficients, of the p F q solutions {/, g} of (80) and their derivatives. 

Based on these observations, the derivation of the solutions (36), (51) and (64) for the BHE, CHE 
and GHE equations (20), (41) and (58), performed in the previous sections, can be reformulated entirely, 
shortcutting the Abel equation step, and resulting in solutions free of integrals as follows. 

The composition of (74) with itself is equal to the identity transformation. 
12 In the easier case co = 0, the equation for p is: p" = cip' + ci'p. 

13 For rewriting equations in normal form, see the Appendix. 

14 An equation is of p F q type if it admits solutions of the form (8); these solutions can be computed systematically [18]. 
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5.1 Solution free of integrals for the BHE (20) 

Applying transformation (74) to the 1F1 equation (33), derived from the AIR (31) in sec. 2, we obtain 

y"= (2(T-x) + -)yi + 2(T + a)xy (81) 



x , 

This equation is already of Heun type, and by rewriting it in normal form 15 , that is, changing further 



y -> exp 



J±- x + rdx^y (82) 



we directly obtain the two parameter BHE (20), derived from the AIR (15). Hence, combining these two 
transformations (74) and (82), we also directly obtain the solution (36) computed in sec. 2 for the BHE (20). 

We note here that the particular form (33) of the 1F1 equation derived in sec. 2, which has only one 
irregular singularity at infinity, has the important feature that under the transformation (74), (33) gains 
one regular singularity at the origin. Hence, the resulting equation (81) is not a 1F1 equation anymore 
but a 2-parameter biconfluent Heun equation. Both the augmentation in the number of singularities under 
transformation (74) and the change in type from 1F1 to BHE do not happen with all 1F1 equations. 

Now, since this BHE equation (81) was obtained from a 1F1 equation using the transformation (74), as 
explained, the derivative p = y' of the solution to (81) also satisfies a 1F1 equation. According to (76), the 
equation for p = y' associated to (81) is 

„ 2{l + rx-x 2 ) , 2{l + TX-(T + a) X z) 

P =— L P 1 t. -p (83) 

x x z 

The solution to this 1F1 equation can be expressed in terms of the Kummer functions M and U [25] as 

p = xc*(™> (M {^-^-, \, (x + .r) 2 ) d + U \, {x + a) 2 ) C 2 ) (84) 

Substituting this solution into (79), we obtain the solution to the BHE (81), and further applying the 
transformation (82), we directly obtain the general solution, free of integrals, for the BHE (20) of sec. 2, as 



+ ((r 2 - a 2 2) M (j-=^ ~ 1, \, (x + ^ — A(x) M (^^, \, (x + a) 2 ) ) C^j (85) 

where A (x) = a 2 + r 2 + 2 (2 x 2 + (3 a — r) x — a t — 1) and a 2 ^ t 2 . Symbolic computation input for verifying 
this solution is found in the Appendix. 

This approach to the solution of the BHE (20) is straightforward, clearly simpler than the calculations 
presented in sec. 2. As shown in the following two subsections, this simpler approach also leads to solutions 
free of integrals for the CHE (41) and the GHE (58) discussed in sec. 3 and sec. 4. We note, however, 
that without the knowledge of the equations of Heun and p F q type being linked, or of the form of the 
transformation relating them 16 , both of which were derived in the previous sections from the connection 
Heun <-> Abel, the existence of the straightforward mechanism used in this section is not evident. 

15 For the formula to write linear equations in normal form, see the Appendix. 

16 In the BHE case, these arc (20), (33) and (36). 
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5.2 Solution free of integrals for the CHE (41) 

Like in the BHE case, a solution to the CHE (41), equivalent to (51) and free of integrals, can be obtained 
by first applying (74) to the 1F1 equation (49), derived from the AIR (47) in sec. 3, leading to a CHE, 

l + 2A( I -l)(l-r-,) 2A»(,-l)(r + g ) 
x (x — 1) x z 

Rewriting this equation in normal form directly results in the CHE (41) discussed in sec. 3. So, the problem 
now is the computation of solutions to the CHE (86). As in the BHE case, we know, by construction, that 
when y satisfies (86), p = y' satisfies a 1F1 equation. According to (76), the equation for p is 

pl , = 3-2A(l + T) + (2A(r + 2)-l)x-2Ax 2 ^ (g?) 

2 A 2 [x - l) 3 (r + a) - 2 {x - 1) (x 2 - 2 x + r + l) A - 1 - x 
+ x 2 {x-lf P 

The general solution to this equation can be written in terms of Whittaker functions M and W [25] as 

(x- iwH^-s/ 2 

p=± J —, (M(u, v, 2 Ax) d+W(n, v,2Xx) C 2 ) (88) 

e x 

where fj, = X(l — <r) + 1/2 and v = A \J t 2 — 2 a + 1. Substituting this solution into (79) leads to a solution 
free of integrals for the CHE (86), from where the solution to its normal form (41) of sec. 3, is 

y= -f=L= ( (W(/i, v, 2\x) + \ (r-a)W(M-l, i/, 2Ax))C 2 
va;- l \ 

+ (\(t + <t)M{ji, is, 2 Xx) + ((1 - a) X - u)M(n - 1, v, 2Ax))Ci^j (89) 

Symbolic computation input for verifying this solution is found in the Appendix. 

5.3 Solution free of integrals for the GHE (58) 

The derivation done in sec. 4 of a solution for the GHE (58) can be reformulated as in the BHE and CHE 
cases. Applying (74) to the 2F1 hypergeometric equation (63), we obtain 

„ _ (2 A + 1) x 2 - 2 ((2 A + 1) a + r) x + 2 a 2 A + (2 r + 1) a , 2 (a - x) A (t + a) y 
V x(x-l)(a-x) V x ^ x -lf m 

which is a GHE with four regular singularities at {0, 1, a, 00}. Rewriting (90) in normal form using 



IP 



, ;2{x-a)T-2{a- X y A + a{2x-l)-x 2 ,\ 

V -> CX P / — " ^7 TY7 v dx y ( 91 ) 

1 2x (x — 1) (x — a) I 



we obtain the 4-paramctcr GHE (58) of sec. 4. The solution (64) presented in sec. 4 for this equation is 
identical to the composition of these two transformations (74) and (91). As in the previous subsections, 
using (76) we compute the equation for p = y' associated to (90), which, as discussed, is by construction 
a 2F1 equation when written in normal form. Solving for p, from (79), we obtain the solution to (90), and 
applying (91) to it, we obtain a solution free of uncomputed integrals for the GHE (58), 
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\/x - a y 

0. f (T - S - A)(T - S + A - """" - a F,(S + A - T + 1,S - A - T + 2, 2 (1 - T); «) 



+ ( T -^) ((« A - T + r ) ^ 1/2_T + (T - E - A) a; 3 / 2 - T ) 2 Fi(E + A - T, E - A - T + 1; 1 - 2T; x)^j 

+ c 2 r (T+s+A)(T+£ - A +i) ^ 



T + ((aA + T + r) x 1 ' 2 ^ - (T + E + A) z 3 / 2 + T ) 2 Fi(E + A + T, E - A + 1 + T; 1 + 2 T; ir) 

where E and T are denned in (66). This solution is also valid when the GHE (58) admits Liouvillian solutions, 
that is, when a 2 = t 2 , although in this case the solutions (59) and (60) are expressed in simpler manner. 
Symbolic computation input for verifying this solution (92) is found in the Appendix. 



6 Comparison with solutions existing in the literature 

A search in the literature didn't show previous references to a link between Heun and Abel equations as the 
one presented in sec. 1, nor a derivation of solutions to the former equations and confluent cases from the 
knowledge of solutions to the latter. It is nonetheless interesting to compare the solutions for Heun equations 
derived through this link Heun <-> Abel and in sec. 5 with the ones previously presented in the literature. 
For practical reasons, the discussion is restricted to three more recent papers, by Ronveax [9], by Ishkhanyan 
and Suominen [10], and by Shanin and Craster [15], which present sufficiently explicit solutions for the GHE 
(1), similar to the non-Liouvillian solutions and the special Liouvillian cases derived in sec. 4.1 and sec. 5.3. 



6.1 Liouvillian solutions 

In [9], the factorization of Heun's General equation (1) into a form 

(L(x)D + M{x)) (L(x)T) + M(x)) y = (93) 

where D = d/dx and {L, M, L, M} are polynomials, is discussed, and six sets of conditions on the Heun 
parameters, such that this type of factorization is possible, are derived, all leading to solutions of the form 

y = x pl (x - If 2 (x - a) p2 (94) 

for some pi. Since these solutions are Liouvillian, they can be computed systematically, e.g., in a symbolic 
computation environment like Maple or Mathematica, where Kovacic's algorithm is implemented. 

The Liouvillian solutions (59) and (60), derived here from the condition a 2 = t 2 related to the canonical 
form of Abel equations, are also of the form (94). Nonetheless, the conditions for the existence of a factor- 
ization of the form (93) obtained in [9] are less general than the conditions for the existence of Liouvillian 
solutions derived here. For example, if in (58) we change variables y — > exp(y), the resulting Heun equation 
will continue having rational coefficients and the condition a 2 = t 2 will continue assuring that the solution 
admits Liouvillian form, computable using Kovacic's algorithm, even when it won't be of the form (94) 
anymore. On the other hand, if we perform the same change of variables in the equations obtained in [9], 
the resulting equations will be out of reach of the factorization there presented, because M{x) in (93) will 
have an exponential factor, while that method applies only to polynomial forms of M(x). 

We note that, by changing variables appropriately, the Liouvillian solutions derived in sec. 3.1 and 4.1 for 
the BHE (20) and CHE (41), can also be transformed into the form (94), which indicates that a factorization 
like the one discussed in [9] exists also for the BHE and CHE. 
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6.2 Non-Liouvillian solutions 



Non-Liouvillian solutions cannot be computed with Kovacic's algorithm, nor is there such a general algorithm 
for computing them. This type of solutions was presented in the previous sections and is discussed in other 
papers in the literature. 

In [10], an approach restricted to the GHE 17 (1) is discussed. Concretely, after some manipulations, 
the equation satisfied by H' , where H is a solution to the Heun equation (1), is presented. This equation 
for H' can be computed using (76), and is a Fuchsian equation with five regular singularites, located at 
{0, 1, oo, a, q/(a /?)}. So, when qj (a (3) is equal to 0, 1 or a, or approaches oo, the equation has four singular- 
ities and hence both H and H' satisfy a Heun equation. This happens when cither q = 0, q = a (3, q = act f3, 
or a [3 — 0. These four cases are presented in [10], and at first sight, the approach could be compared with 
the one presented here, in sec. 5, where Heun and related confluent equations with the property that H' is 
of the form (8), involving p F q functions, were derived. 

The main difference between the presentation in sec. 5 and that in [10] is that, in sec. 5, the non-local 
transformation (74) directly leads to non-trivial Heun equations in H, such that H 1 can be computed 
systematically 1 because it is of the form (8), and from there we can systematically compute H, using (79). 

On the other hand, in [10] , the condition that H' satisfies a Heun equation of one parameter less is not 
sufficient to compute its value. So, to obtain solutions using this approach, the authors introduce additional 
restrictions on the values of the Heun parameters so that H' admits solutions expressible using 2F1 functions. 
Neither the origin of these ad-hoc restrictions nor a systematic manner of computing them is shown. 

For the first case, q = 0, the additional restrictions suggested in [10] are e = — 1 and q' = aa'/3', where 
{q 1 , a\ [3'} are some functions of the Heun parameters {a, (3, 7, 5, e, q}, together leading to 



This is a case depending on only three parameters {a 1 , (3' , 7}, which happens to be a particular case of the 
4-parameter GHE solved here, in sec. 5.3. That can be seen by using (76) to construct the equation for 
p= H' associated to Heun equation (90) of sec. 5.3, and by noting that its solutions arc of the form (8) 



for some pf, that is, they depend on four parameters {a,a,/3, 7}, not three. 

The other case explicitly discussed in [10] is q — a (3, so a 5-parameter GHE for which a solution is shown 
in terms of an infinite series of Appel functions 18 [28] which, in one case, is shown to be expressible as an 
infinite sum of 2F1 functions, presented in [10] with number (36). It is not evident how to compare this 
formal infinite series solution with the finite-number-of-terms solutions presented here, in sec. 5.3, nor is it 
evident how to impose an additional restriction to the 5-parameter GHE treated in [10] such that the formal 
infinite series expansion terminates. Solutions to the GHE in terms of infinite series of 2F1 functions are also 
known in the literature [1], although the authors of [10] make the point that the 2F1 functions entering the 
formal series they present have a different behavior than those shown in [1]. 

An approach based on removing "false" singularities 

A referee has also pointed to [15], a very interesting paper by Shanin and Craster, published in 2002 but 
actually submitted during 2000, where a more thorough approach to solving GHE and confluent CHE 
equations (actually, the approach is for all linear equations with "false" singularities), is presented. 

The first idea in [15] consists of determining relations (constraints) between the Heun parameters such 
that one of the regular singularities is "false". That will happen when the corresponding characteristic 
exponents differ by an integer but, also, no logarithmic term appears in the local expansion of the solution. 
Although such an approach is entirely different from the one developed here, where the solvable Heun 
equations are derived from the single Abel AIR (11), it is remarkable that the BHE (20), CHE (41) and 
GHE (58) resulting from this link Heun <-> Abel do have one such false singularity. 

17 The BHE, CHE or other confluent cases are not discussed in [10]. 

18 The Appel function is a formal extension of the 2F1 function to two variables, expressed as a double infinite series. 



H' = x 2 F 1 (c/,/3'; 7 + 2;x) 



(95) 



p = H' = x Pl (x-l) P2 (x- a) P3 2 Fi (a, (3; 7; x) 



(96) 
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In [15], there is no explicit discussion of the BHE (4), but the approach seems applicable to that case 
too, and, as is the case here, the approach does not seem to be applicable to the DHE (5) or the THE (6). 
The derivation of GHE and CHE solvable equations in [15] is systematic but not as straightforward as the 
one-step derivation shown here in sec. 5.2 and sec. 5.3, exploring non-local transformations. Concretely, in 
[15], determining the value of the accessory parameter q in (1), such that the equation has a false singularity, 
requires using p F q identities and solving recurrence relations for the coefficients of series expansions, to assure 
there is no logarithmic term in the solution. 

The second idea presented in [15] is that, when the Heun equation has a false singularity, its solution can 
be expressed as a linear combination (with finite number of terms and constant coefficients) of p F q functions, 
determined by exploring isomonodromy mappings. Although finding the appropriate isomonodromy is a 
powerful idea, and the approach is systematic, quoting the authors of [15] (p. 628): "The procedure of finding 
an appropriate isomonodromy mapping described is quite complicated. In the relatively simple examples that 
we construct we can pose an ansatz for the form of the mapping to within several unknown parameters; 
these are found by direct substitution. ". The actual procedure to determine these unknown parameters is 
systematic, but not so straightforward; quoting the authors (p. 630): "The simplest way to do this is to 
substitute the linear combination, say U + CV directly into the Heun equation and then use known recursion 
formulae for hypergeometric functions and their derivatives. Tedious calculations show that ... ". 

In contrast, the approach presented here, in sec. 5, directly leads, by construction, to the three multipa- 
rameter BHE (81), CHE (86) and GHE (90) equations having false singularities as well as to the exact form 
of the linear combinations (with non-constant coefficients) of p F q functions that solve these equations. 

Furthermore, here, in sec. 5, the use of non-local transformations links p F q equations also to other linear 
equations, with more singularities than those of the Heun families, where again the linear combinations (79) 
solving all these equations involve non-constant coefficients. Although these cases can in principle be treated 
by finding isomonodromies, that approach may result non-practical as soon as the number of singularities 
or the number of parameters involved increases. For example, in perhaps the simplest case, departing from 



qFjX ; a; x) 

X — K 

the p F q equation satisfied by this expression is 



(97) 



i a 2 \ , (x — k — a) y 

y"=[ — + )y+—, 9 8 

X K— X ) X\X— K) 



Applying now the non-local transformation (74) we obtain 



, ^ — 1 1 1 \ , (x — k — a) y , . 

y" = + + y' + (99) 



x (x — k) 

This equation has three regular singularities at {0, k, a+ n} and one irregular singularity at oo; they are 
all irreducible, and therefore (99) does not fit into any of the five Heun classes represented by (2-6) (it 
belongs to an "upper" class). According to sec. 5, by construction, the equation satisfied by p = y' admits 
systematically computable [18] solutions of the form (8) which, when plugged into (79), lead to the following 
solution to (99): 

y = d x a {a „Fi( ;a;x)-(x-K) F 1 (;a+ 1; a:)) (100) 
+ C* 2 ((a - 2) ((1 -a)n + ax) F 1 ( ; 2 - a; x) + x (x - k) Fi( ; 3 - a; x)) 

This approach, as described in sec. 5, works just as straightforwardly when we start with p F q functions 
more general than (97), while through that process the equation resulting from applying (74) can be made 
to depend on more parameters and have more singularities. Even so, by construction, the exact linear 
combination (with non-constant coefficients) of p F q functions solving the resulting equation is always given 
by (79). Contrasting with that, depending on the starting p F q expression to be used in place of (97), the 
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construction of the same solvable cases and computation of their solutions using the approach presented in 
[15] can be really complicated. 

On the other hand, an important generalization presented in [15] is that it provides a recipe for computing 
the isomonodromies and constructing the related Heun equations having as solutions linear combinations of 
p F q functions involving more than two terms. 



7 Discussion 

In sec. 2, 3 and 4, solutions in terms of p F q functions were derived for families of the Heun equations BHE, 
CHE and GHE. The approach links linear equations with four regular singularites (and related confluent 
cases) to linear equations with three regular singularities (and related confluent cases), by linking both types 
of linear equations to the canonical forms of the Abel AIR class of non-linear first order equations. The link 
AIR <-> p F q is developed in [17], and the link AIR <-» Heun is presented in sec. 1. This link also provided a 
natural way to determine the special Liouvillian cases of the Heun solutions to the BHE, CHE and GHE here 
treated, and permits studying Abel equation problems by reformulating them in terms of linear equations. 

In sec. 5, that approach is shown to be equivalent to performing the non-local transformation (74) on 
the p F q equations (33), (49) and (63), leading to Heun equations with two important properties: 1) further 
auxiliary equations which can be derived from them for p = y' are of p F q type; 2) when written in normal 
form, these Heun equations obtained using (74) are identical to the BHE (20), CHE (41) and GHE (58) 
solved in the sections previous to sec. 5. This approach leads in a simpler manner to the same solutions (36), 
(51) and (64), and also to the equivalent forms of these solutions free of integrals, (85), (89) and (92). All 
the solutions presented were verified for correctness using symbolic computation software. 

Besides the presentation in sec. 2, 3 and 4, the existence of a connection between the "Heun and related 
confluent equations" on the one hand and the "AIR (11) and the different possible multiplicities of its roots 
Pi" on the other hand, can be seen more straightforwardly by transforming not (14) but (11) into a linear 
equation 19 , resulting in 

2/ »=p- + -^ + -^ + -^V + (s °" + ro)( T x) — 2 v (ioi) 

\x-a x-p 3 x-p 2 x-pij ( x - p 3 ) {x - p 2 ) {x - Pi) 
This is a Heun equation with its four regular singularites at {pi, p 2 , p 3l a}, where Ri + R 2 + R3 = —3, 

D P2 2 - (Sl +P3 + Pl)p2 + P3pl - n 

R2 = 7 77 s (1U2J 

(Pi - Pi) (pi - Pz) 

and {Ri,Rz} are obtained from R 2 multiplying by —1 and respectively swapping p 2 <-> p\ and p 2 <-> p 3 . 
Through the confluence processes which coalesce singularities in (1), generating the CHE and BHE confluent 
equations, one coalesces the corresponding singularities pi of (101), generating the same type of confluent 
equations, and that is equivalent to having multiple roots pi in the AIR (11). 

By rewriting (101) in normal form (see (58)), the number of irreducible parameters involved is shown to 
be four instead of six as in (2). That explains the restrictions on the Heun parameters of the BHE, CHE 
and GHE families discussed in the previous sections. In the three cases, one parameter is fixed and another 
one is dependent on those remaining. 

Different from the BHE, CHE and GHE cases, in the case of the DHE (5) and THE (6) the approach 
considered in this paper does not lead to new solutions. That can be seen by applying to (101) the DHE 
and THE confluence processes [2], in both cases arriving at equations already of p F q type. That status of 
things is somewhat expected: the AIR class is generated from the three canonical forms (14) and these are, 
in their general form, already linked to GHE, CHE and BHE families. 

Independent of the possibility, developed here, of expressing the solutions to the BHE (20), CHE (41) 
and GHE (58) normal forms without introducing "Heun functions", these functions have been developed 
consistenly during the last years and will most certainly form part of the standard mathematical language 

19 For that purpose, apply first {x <-> y} to (11) at {s2 = l,r-2 = —a}, then apply (12) to the resulting Riccati equation. 
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in the near future. That can be inferred from the relevance of Heun equations in applications. In this 
framework, the results of this paper could be seen as the identification of multi-parameter special cases of 
Heun functions of the BHE, CHE and GHE types, respectively admitting the integral representations (36), 
(51) and (64), and the linear combinations of p F q functions with variable coefficients (85), (89) and (92). The 
mathematical properties and the relevance of these special cases in applications require further investigation. 

This link between Heun and p F q second order linear equations through Abel non-linear equations of first 
order seems to be the simplest case of a link between linear equations with N and N-l singularities, through 
"Abel AIR like" equations, for which the numerator of the right-hand-side has degree N-l at most 20 . For 
example, if instead of (11) we depart from 

v' = (y-p^iy-piUy-P^y-P,) (103) 

(s 2 x z + si x + s j y + r 2 x z + n x + r 

that is, an equation with structure similar to the AIR (11) but whose numerator of the right-hand-side is 
of degree four, then by applying {x <-> y} to obtain a Riccati equation and transforming the latter into a 
second order linear equation, we obtain an equation similar to (101) but with five regular singular points, 



y » =(J—+ Ri + R * + R2 + Rl )y> + ( SO x + r Q )(a-x) 

\x-a x-p 4 x-p 3 x-p 2 x-px) (x - pif (x - p 3 f (x - p 2 f (x - Pl f 

which, together with its confluent cases, can be linked through (103), this time to the Heun equations, using 
the same approach presented in the previous sections relating Heun to p F q equations. 

Analogously, n th order (n > 2) linear equations in y(x) can also be reduced to "Riccati like" non-linear 
equations of order n — 1, due to their invariance under scalings of y. It is therefore reasonable to expect 
that a link equivalent to the one discussed in this work also exists between linear equations with N and N-l 
singularities in the n th (n > 2) order case. 



Appendix 

In sec. 2, 3 and 4, solutions were derived for the equations in normal form BHE (20), CHE (41) and GHE 
(58). Given a second order linear ODE 

y" + c lV ' + c y = Q (105) 
where the Cj = Cj (x) , the corresponding normal form, 

y"+ (c -( Cl 2 + 2 Cl ')/4)?/ = (106) 

is obtained by changing y — > cxp(— fa dx/2) y. 

Regarding Heun equations, one advantage of the normal form is that the general or confluent type of the 
equation is evident in the partial fraction decomposition of the coefficient of y (see eqs. (2) to (6)). Also, 
two different equations related by y — > P(x) y have the same normal form and recognizing this equivalence is 
relevant for computational purposes. On the other hand, for different reasons, special functions are frequently 
defined as solutions to equations in canonical form. This appendix relates the normal and canonical forms 
of the BHE, CHE and GHE, expressed in terms of irreducible parameters 21 , following the notation of [2], 
thus permitting a simple translation of the results presented. 

20 For Abel equations of the first kind, N = 3; for Abel equations of the second kind like (9), N < 3. 
21 A few of the equations shown in the paper are repeated here for ease of reading. 
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The Biconfluent Heun equation 

The BHE canonical form is given in terms of four constant parameters {a, (3, 7, 6} by 

y'' + (^-P-2 X )y' + {,-a-2- S+ ^ + x a ^ )y = 0-, (107) 
The BHE in normal form (4) in terms of four parameters {B,C,D, E} is 

y"+ ( -x 2 + Bx + C+- + ^r )y = (108) 



The BHE normal form (20) solved in sec. 2, there written in terms of two parameters {cr, r}, is 

y»-( x 2 + 2ax + T 2 + ^- + ^y = (109) 

The parameters {B, C, D, E} in (108) are related to {a, r} by 

B = -2a, C=-D 2 , D = -t, £=-3/4 (110) 
The parameters {a, (3, 7, 5} in (107) are related to {B, C, D, E} by 

a 2 = -AE+l, 0=-B, ~i = B 2 /A + C, 5 = -2D (111) 
So the parameters {a, (3, 7, 6} are related to {cr, t} by 

a 2 =4, (3 = 2a, 7 = cr 2 - t 2 , 5 = 2t (112) 

At these values of {a, f3,j,5}, for a = ±r, (107) admits Liovillian solutions and for a 2 ^ r 2 the solution is 
obtained from (85). 

The Confluent Heun equation 

The CHE canonical form is given in terms of five constant parameters {a, (3, 7, 5, r]} by 

\ x x — 1 / 2x (x — 1) 

The CHE in normal form (3) in terms of five parameters {A, B, C, D, E} is 

»"+(- 4 + ! + ; £ t + # + <^if)» = (114) 

The CHE normal form (41) solved in sec. 3, there written in terms of three parameters {A, cr, r}, is 



„ , ., 2 (a- 1)A 2 -rA + 1/2 tX-1/2 (t 2 - 2 a + 1) A 2 - 1/4 3 , 

^"-^ 2 + - — x L + -^~ + - ^ - + I{x-Tf ]y = (115) 

The parameters {A, B, C, D, E} in (114) are related to {A, a, r} by 



,4 = -A 2 , B = 2(1 -cr) A 2 +r A- 1 C*=i-rA, £> = 1 + (2 cr - r 2 - l) A 2 , £ = -3/4, (116) 
The parameters {a, (3, 7, <5, 77} in (113) are related to {A, B, C, D, E} by 
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a 2 = -4 A, (3 2 = -AD+l, 7 2 =4 7 -4£-3 5 = C + B - a, rf = ~ - B - (3 (117) 
So the relation between {a, (3, 7, 5, 77} and {A, a, r} is 

a 2 =4A 2 , /3 2 = 4 (l — 2 (T + r 2 ) A 2 7 2 = 4 7 , (5 = 2 (1 - cr) A 2 - a, 77 = 2 (a - 1) A 2 - r A - /3 (118) 

At these values of {a, (3, 7, 6, rj}, the CHE (113) admits Liouvillian solutions for a = ±t, and for a 2 ^ t 2 
the solution is obtained from (89). 

The General Heun equation 

The GHE canonical form is written in terms of seven constant parameters {a, (3, 7, S, e, a, q} as 

„ / 7 S e \ , a 13 x — q , 

y"+[ 1 + — 7 + h/ + ? \ w r y = 119 

\x x—1 x — a) x(x — l)(x — a) 

where j + 5 + e = a + (3+l and a^0,a^l. In [2], the numerator of the coefficient of y of this equation 
is written as a (3 (x — h). The notation in (119) is the one used in [1], was apparently first adopted in [27], 
and has the advantage that one can take a = (or (3 — 0) without eliminating the term in y completely 22 . 
The GHE in normal form (2) in terms of six parameters {a, A, B, D, E, F} is 



1 A B A + B D E F ■, 

y +[- + — 7 + — + 7 — 772 + 7 72 )v = 120 

x x-l x-a x 2 (x-1) (x-a)J 
The GHE normal form solved in sec. 4 (see (58)), written in terms of four parameters {a, cr, t, A}, is 



„ _ (2a 2 (a- 1) A 2 - 2aa(2a- 1) A+ (2t 2 - l/2)a + r + 1/2 

y — (121) 

\ ax 

2 (a (a- l) 2 A 2 - o (2 a - 1) (a - 1) A + (r - 1/2) ((t + 1/2) a - t)) r - a + 1/2 



(a — 1) (x — 1) a (a — 1) (2 — a) 

Q 2 A 2 -2acrA + T 2 -l/4 (a - l) 2 A 2 - 2 g (a - 1) A + t 2 - 1/4 3 \ 

+ ~ 2 + (x~lf + A{x-a) 2 J V 

The parameters {A, B, D, E, F} in (120) are related to {a, cr, t, A} by 

A = -2a(a - 1)A 2 + 2(2q - l)gA-2r 2 - T + 1/2 + ^, 

CI A 



T ■ 



a/2 



B = 2a(a-l)A 2 -2(2a-l)crA + 2T 2 + j . 

D = -a 2 A 2 + 2acrA-T 2 + l/4, (122) 
E = -(a-l) 2 A 2 + 2(a-l)crA-T 2 + l/4, 
F = -3/4 

The parameters {a, (3, 7, S, e, c/} in (119) are related to {a, A, B, D, E, F} by 

22 The differences in notation with respect to [2] are the coefficient of y and the use of {7, <5, e, a, (3, q} in place of {a, (3, 7, 8, rj, h}. 
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7 2 = -4L> + 2 7 , 5 2 = -4E + 2S, e 2 = -AF + 2e, 
a 2 = ((S + E + 7-l)a- (7+ 2 ^ -^ + ( tt -l)g + Bj 4 

(3 = j + 5 + e-a-l, q =^±±L- a A 
So the parameters {a, (3, 7, 5, e, q} are related to {a, a, t, A} by 



7 2 


= 2(2a 2 A 2 -4acrA + 2r 2 + 7 )-l 






= 2 (2 (a - l) 2 A 2 + 4 a (1 - a) A + c> + 2 t 2 ) - 1 




6 2 


= 2e + 3 




a 2 


= 2 a (1 - a) A 2 + 2 (2 a - 1) cr A + (7 + (5 + e - 1) a - 


2T 2 -(75+l + (7 + (5)e)/2 


a 


= 7 + 5 + e — a — 1 




Q 


a 5 + 4t 2 — l) 

- 2aA(a(a 1) A cr(2a 1))+ w 7 


76 1 

+ T + T+ 2 



At these values of {a, j3, 7, S, e, q}, for cr = ±r the GHE (119) admits Liouvillian solutions, and for cr 2 7^ r 2 
the solution is obtained from (92). 

Verifying solutions using symbolic computation 

In presentations like this one, where equations and solutions involving many parameters and non trivial 
special functions are involved, it is of use to be able to verify the correctness of the solutions derived, in 
some way alternative to the one presented. For that purpose, the input lines, written in the Maple symbolic 
computation syntax, for the BHE (20), the CHE (41) and the GHE (58) equations and their respective 
solutions (85), (89) and (92), are given, so that they can be copied from the online version of this paper. 
The 2-parameter BHE (20) is written in Maple syntax as 

> BHE := diff (y(x) ,x,x) - (x~2 + 2*sigma*x + tau~2 + tau/x + 3/4/x~2) *y (x) = 0; 

and its solution (85) is written as 

> BHE_sol := y = exp(-sigma*x-l/2*x~2)/x~(l/2)/(x+sigma)*((Lambda*KummerU(l/4*tau~2 

> - l/4*sigma~2 , 1/2 , (x+sigma) ~2) -4*KummerU(l/4*tau"2-l/4*sigma~2-l , 1/2 , (x+sigma) ~2) ) 

> * _Cl+( (tau~2-sigma~2-2) *KummerM(l/4*tau~2-l/4*sigma~2-l , 1/2 , (x+sigma) ~2) -Lambda 

> * KunimerM(l/4*tau~2-l/4*sigma~2,l/2, (x+sigma) ~2) ) *_C2) ; 

> Lambda := sigma~2+tau~2+4*x~2+2*(3*sigma-tau)*x-2*sigma*tau-2; 

After entering these lines in a Maple session, to verify this solution one can use the Maple odetest command, 
as in > odetest ( BHE_sol, BHE ) ; which returns zero, confirming that the solution cancels the equation. 
The 3-parameter CHE (41) is written in Maple syntax as 

> CHE := diff (y(x) ,x,x) = ((-l+2*tau*lambda)/(2*x-2)+l/2*(l+(-4+4*sigma) 

> * lambda~2-2*tau*lambda) /x+1/4* (-l+(4-8*sigma+4*tau~2) *lambda~2) /x"2 

> + 3/4/(x-l)~2+lambda~2)*y(x) ; 

and its solution (89) is written as 

> CHE_sol := y(x) = 1/ (x-1) " (1/2) * ( ( (tau+sigma) *lambda*WhittakerM(mu,nu,2*lambda*x) 

> + ( (-sigma+1) *lambda-mi) *WhittakerM(-l+mu,mi,2*lambda*x) ) *_C1+ (lambda* (tau-sigma) 

> * WhittakerW(-l+mu,nu,2*lambda*x)+WhittakerW(mu,nu,2*lambda*x))*_C2) ; 

> mu := l/2-lambda*sigma+lambda; nu := lambda* (-2*sigma+l+tau~2) " (1/2) ; 
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The 4-parameter GHE (58) is written in Maple syntax as 

> GHE := diff (y(x) ,x,x) = ( (2*a~2* (a-1) *Delta~2-2*sigma*a* (2*a-l) *Delta 

> + (2*tau~2-l/2) *a+tau+l/2) /x/a-2* (a* (a-1) ~2*Delta~2-sigma* (2*a-l) * (a-1) *Delta 

> + (tau-l/2)*((tau+l/2)*a-tau))/(x-l)/(a-l)+(tau-a+l/2)/a/(a-l)/(x-a) 

> + (Delta~2*a~2-2*a*sigma*Delta+tau~2-l/4) /x~2+ ( (a-1) ~2*Delta~2-2* (a-1) *sigma 

> * Delta+tau~2-l/4) /(x-1) "2+3/4/ (x-a) ~2) *y(x) ; 

and its solution (92) is written as 

> GHE_sol := y(x) = 1/ (x-a) ~ (1/2) * (x-1) ~ (1/2+Sigma) * ( (1/2* (Tau-Sigma-Delta) 

> * (-Sigma+Delta-l+Tau) * (x~ (5/2-Tau) -x~ (3/2-Tau) ) *hypergeom( [Sigma-Delta+2 

> - Tau, Sigma+Delta-Tau+1] , [2-2*Tau] ,x)+hypergeom( [Sigma+Delta-Tau, Sigma-Delta 

> + 1-Tau] , [l-2*Tau] ,x) *( (Tau-Sigma-Delta) *x~ (3/2-Tau) +x~ (1/2-Tau) 

> * (-Tau+a*Delta+tau) ) * (-1/2+Tau) ) *_C1+_C2* (1/2* (Sigma+Delta+Tau) * (Sigma-Delta 

> + l+Tau)*(-x~(5/2+Tau)+x~(3/2+Tau))*hypergeom( [Sigma-Delta+2+Tau, Sigma 

> + Delta+Tau+1] , [2+2*Tau] ,x)+( 1/2+Tau) *( (-Tau-Sigma-Delta) *x~(3/2+Tau) 

> + x~( 1/2+Tau) *(Tau+a*Delta+tau)) 

> * hypergeom( [Sigma+Delta+Tau, Sigma-Delta+l+Tau] , [l+2*Tau] ,x) ) ) ; 

> Sigma := sqrt ( (a-1) ~2*Delta~2-2* (a-1) *sigma*Delta+tau~2) ; 

> Tau := sqrt(a~2*Delta~2-2*a*sigma*Delta+tau~2) ; 
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